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ABSTRACT 

Context. Interpretation of solar polarization spectra accounting for partial or complete frequency redistribution requires data on various col¬ 
lisional processes. Data for depolarization and polarization transfer are needed but often missing, while data for collisional broadening are 
usually more readily available. Recent work by Sahal-Brechot and Bommier concluded that despite underlying similarities in the physics of 
collisional broadening and depolarization processes, relationships between them are not possible to derive purely analytically. 

Aims. We aim to derive accurate numerical relationships between the collisional broadening rates and the collisional depolarization and polar¬ 
ization transfer rates due to hydrogen atom collisions. Such relationships would enable accurate and efficient estimation of collisional data for 
solar applications. 

Methods. Using earlier results for broadening and depolarization processes based on general (i.e. not specific to a given atom), semi-classical 
calculations employing interaction potentials from perturbation theory, genetic programming (GP) has been used to fit the available data and 
generate analytical functions describing the relationships between them. The predicted relationships from the GP-based model are compared 
with the original data to estimate the accuracy of the method. 

Results. Strongly non-linear relationships between the collisional broadening rates and the depolarization and polarization transfer rates are 
obtained, and are shown to reproduce the original data with accuracy of order 5%. Our results allow the determination of the depolarization 
and polarization transfer rates for hyperfine or fine-structure levels of simple and complex atoms. 

Conclusions. In this work we have shown that by using a numerical approach, useful relationships with sufficient accuracy for applications are 
possible. 
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1. Introduction 

Interpretation of polarization in solar spectral lines (the so- 
called second solar spectrum) accounting for frequency redis¬ 
tribution requires understanding of collisional processes occur¬ 
ring in the solar atmosphere. In particular, data for depolariza¬ 
tion and polarization transfer by collisions with the most abun¬ 
dant perturber, hydrogen atoms, are needed, but often missing 
and thus neglected in modelling. On the other hand, data for 
collisional broadening due to hydrogen collisions are usually 
more readily available. Thus, possible relationships between 
these processes to enable accurate and efficient estimation of 
depolarization data would be useful in understanding the for¬ 
mation of the second solar spectrum. 

During the 1990s, Anstee, Barklem and O’Mara (ABO) de¬ 
veloped a semi-classical theory for collisional line broadening 
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by neutral hydrogen (Anstee 1992; Anstee & O’Mara 1991, 
1995; Anstee et al. 1997; Barklem 1998; Barklem & O’Mara 
1997; Barklem et al. 1998). That theory was later extended by 
Derouich, Sahal-Brechot and Barklem (DSB) to the depolar¬ 
ization and polarization transfer by collisions with neutral hy¬ 
drogen (see for example Derouich et al. 2003a and Derouich 
2004). The same potentials, based on Rayleigh-Schrodinger 
perturbation theory in the Unsold approximation, so called 
RSU potentials, were used in both cases. In addition, the same 
time-dependent Schrodinger equation was solved to obtain the 
scattering S-matrix; an essential ingredient in the determination 
of the both collisional depolarization and broadening rates. 

A great advantage of the ABO and DSB methods is that 
they are general, i.e. not specific to a given perturbed atom/ion, 
and therefore adaptable to any neutral or singly ionised atom. 
The general, semi-classical methods of ABO and DSB give 
results in good agreement with accurate but time consuming 
quantum chemistry calculations to better than 20% for solar 




2 


Derouich et al.: Collisional effects 


temperatures (T ~ 5000 K). For instance, at T =5000K, the dif¬ 
ference between results for depolarization rates from the gen¬ 
eral semi-classical theory (Derouich et al. 2003a) and quantal 
results is 11% for Mgi 3p 1 / J i, 8% for Cai 4p l P\, and 13% 
for Na I 3p 2 / J 3/2 • In addition, Derouich et al. (2006) showed 
in their Fig. 6 that the error bar on their collisional rates deter¬ 
mination is well located within the expected error bar on the 
value of the solar magnetic field. They showed that the error 
bar for the magnetic field determination due to the typical un¬ 
certainty in DSB collisional depolarization rates is similar to 
the error bar for the magnetic field value due to a polarimetric 
uncertainty of ~ 10 4 , which corresponds to the uncertainty at¬ 
tainable with polarimeters such as ZIMPOL (Zurich IMaging 
POLarimeter) and THEMIS (Heliographic Telescope for the 
Study of the Magnetism and Instabilities on the Sun). 

The underlying similarity between collisional depolariza¬ 
tion and broadening processes as calculated by DSB and ABO 
suggests the possibility of establishing relationships between 
them. The possibility of obtaining such relationships was ex¬ 
amined from a theoretical viewpoint by Sahal-Brechot and 
Bommier (2014). However, they concluded that there is no 
possible purely analytical relationship between the collisional 
depolarization rate and the collisional broadening coefficient, 
even if the same basic methods and computer codes can be used 
for obtaining the interaction potentials and collisional scatter¬ 
ing matrix. In this work we aim to provide such relationships 
between the results of the ABO and DSB methods through a 
numerical approach. 

The paper is organized as follows. Section 2 gives a brief 
review of the basic definitions and notations, and describes the 
collisional data employed in the context of this work. Section 
3 explains the optimization approach used to obtain the unified 
numerical model. The results for simple atoms without hyper- 
fine structure, simple atoms with hyperfine structure and com¬ 
plex atoms are presented in section 4. Finally, the conclusion 
of the paper is presented. 

2. Background and theory 

The theory of partial or complete redistribution of polarized ra¬ 
diation in a magnetized plasma, which permits to determine the 
Stokes parameters, must account for the effects of collisions 
(e.g. Stenflo, 1994; Bommier, 1997a,b; Landi Degl’Innocenti 
et al., 1997, Casini et al. 2014). The resultant effect of the col¬ 
lisional depolarization and the collisional broadening by neu¬ 
tral hydrogen is accompanied by a variety of interesting spec- 
tropolarimetry effects. Collisions act mainly with three distinct, 
yet correlated, effects; depolarization and polarization trans¬ 
fer, line broadening, and frequency redistribution in scatter¬ 
ing events (partial destruction of the correlation between in¬ 
cident and scattered frequencies). The theory of collisions and 
their relationship to the theory of polarization in spectral lines 
is expounded in detail elsewhere (e.g. Stenflo 1994, Landi 
Degl’Innocenti & Landolfi 2004). Here we simply recall im¬ 
portant relations for our work and define our notations. 

ABO calculations are performed in the standard dyadic ba¬ 
sis (Anstee 1992; Anstee & O’Mara 1991, 1995; Anstee et al. 
1997; Barklem 1998; Barklem & O’Mara 1997; Barklem et al. 


1998). Let us consider the broadening of the line profile be¬ 
tween the initial electronic state | nl\) and the final state \nl 2 ) 
where n is the principal quantum number and /1 and l 2 are the 
orbital angular momentum quantum numbers of the states |n/i) 
and ln/ 2 ), respectively. We denote by ////, and /«/, the projec¬ 
tions of the orbital angular momentums l\ and l 2 on to the in¬ 
teratomic axis, which is taken as the axis of quantization (see 
ABO papers for more details). Fig. 1 shows examples of the 
collisions resulting in the broadening of the line profile for the 
transition between the initial electronic state |n/i) and the fi¬ 
nal state In/ 2 ). Note, the ABO theory neglects the effects of the 
spin, and thus fine structure is ignored. 
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Fig. 1. As in the ABO theory, we consider two electronic states 
|n/[) and \nl 2 ) connected by a radiative transition. The colli¬ 
sional broadening of the line pofile between the levels \nl\) and 
In/ 2 ) is a linear combination of all the collisional rates between 
the Zeeman sublevels inside the state |n/i) and the collisional 
rates between the Zeeman sublevels inside the state In/ 2 ). Note 
that the level spacings are not to scale. 


In order to calculate the depolarization rates of the fine 
structure level \nlJ), DSB used the basis of irreducible tenso- 
rial operators and the atomic states are described by the density 
matrix elements p k q (J) where the tensorial order 0 < k < 2J 
and the coherence -k < q < k. More details about the 
physical meaning of the tensorial order k and the coherence 
q can be found, for instance, in Omont (1977), Sahal-Brechot 
(1977), and Landi DeglTnnocenti & Landolfi (2004). We de¬ 
note by D k (J) the depolarization rates and by <r k the depolar¬ 
ization cross sections. The depolarization rate of the fine struc¬ 
ture level | nlJ) is a linear combination of the purely elastic col¬ 
lisions which occur between Zeeman sublevels | nlJMj) (see 
Fig. 2). Mj is the projection of the total angular momentum J 
on the interatomic axis which is taken as the axis of quantiza¬ 
tion (see DSB papers for more details). 


















Derouich et al.: Collisional effects 


3 



M 

Fine structure level J / 

\ J 

\ 

* nr 


J 


Electronic state / nl > 


where DSB found a similar power law relationship with veloc¬ 
ity: 

= J' and J*f) = ff*_ y ,(vo) W ^ , (3) 

with Aj }J , ~ 0.25. The DSB cross sections, as the ABO cross 
sections, are tabulated for a reference velocity of vo = 10 4 m/s. 

The probability that a collision does not destroy the polar¬ 
ization associated with a 2/t-multipole during a scattering event 
is (e.g. equation 10.19 of Stenflo 1994): 


Fig. 2. The depolarization rate is a linear combination of the 
collisional rates between the Zeeman sublevels inside the state 
| nlJ). Note that the level spacings are not to scale. 


The polarization transfer rates from the \nlJ\) to the level 
\nIJ 2 ) are denoted by C k (J 1 —> Ji)\ cr^ :iJ are the polarization 
transfer cross sections. As it is showed in the Fig. 3, transfer 
of polarization between the levels \nlJ\) and \nlJ 2 ) is a lin¬ 
ear combination of the inelastic collisional rates between the 
Zeeman sublevels of the level \nlJ\) and the Zeeman sublevels 
of the level \nlJ\). Note that, the inealstic collisions which con¬ 
tribute in the polarization transfer rates occur inside one elec¬ 
tronic state. 
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where y E is the elastic collision rate, y N the radiative broad¬ 
ening rate, and D k is the rate of depolarizing collisions. The 
probability that the tensorial order k will be collisionally de¬ 
stroyed during the scattering process is (equation 10.21 of 
Stenflo 1994): 


P 2 


D k 

7n + D k ' 


(5) 


The ratio of these probabilities is 


Pi_ _ 7 e - D k _ 1 ~ r 
Pi D k r 


(6) 


where r = D k /y E . Note, r should not be greater than unity since 
the rate of elastic collisions leading to depolarization may not 
exceed the total elastic collision rate. 

If one considers an isolated spectral line where inelastic 
collisions are assumed negligible, and lower state interactions 
ignored, in the impact approximation the full width at half 
maximum (FWHM) of the collisionally broadened line, y E , is 
equal to the total elastic collision rate y E ; see, e.g., pg 517 of 
Baranger (1962). In line broadening calculations, including the 
ABO theory, it is usual to define the broadening cross section 
such that it gives the half width (HWHM), here denoted w (see 
below). Thus, in such a case, r = D k /y E ~ D k /y E = D k /(2w), 
and so in this work we find it convenient to define: 


R = 


(7) 


where 


r w = D k /w = 2 r. 


( 8 ) 


Fig. 3. Schematic representation of the collisional transfer of 
polarization between different fine structure levels. Note that 
the level spacings are not to scale. 


The depolarization rates and the polarization transfer rates 
are : 


r»oo 

D k (J) = n H v/(v, T)cr k j(v), 

Jo 

(1) 

s>oo 

C k (J — > /') = n H V /(v, T^jAy), 

Jo 

(2) 


The collisional line width is defined in terms of the broad¬ 
ening cross section cr w , by: 

/-*oo 

w = yd 2 = n H v /(v, T)cr w (v), (9) 

Jo 

where /(v, T) is the Maxwellian velocity distribution for a local 
temperature T and the perturber density. ABO found that the 
cross sections have approximately a power law behaviour with 
velocity, such that they define: 

cr H ,(y) = £T w (v 0 ) — 

\> ; 0 


(10) 
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where vo = 10 4 m/s is the reference velocity for which cross 
sections are tabulated and A is the so-called velocity exponent. 
ABO found that A has a limited range of variation about 0.25 
and, as a result, that w has typical resulting temperature depen¬ 
dence of r 038 . 


We note that several works have attempted to make simple 
estimates of the relationship between broadening and depolar¬ 
ization. Smitha et al. (2014) calculated for the case where the 
interaction between particles can be described by a single ten¬ 
sor operator of given rank, using the equations given by Landi 
Degl’Innocenti & Landolfi (2004). They calculated for their 
specific case of a Sc ii line with upper state having J - 2 for 
rank 1 and 2 operators, and found that r = D 1 Ije = 0.1 and 
0.243 (or r„ = 0.2 and 0.486), respectively. Interactions of rank 
1 correspond to the dipolar interaction, and rank 2 correspond 
to dipole-dipole interactions; neither is strictly applicable to 
atom-atom or atom-ion interactions. Faurobert et al (1995) find 
r = D 2 /jc- 0-5 (or r w = 1), for the Sri resonance line. 


The ABO and DSB calculations consider a simple neutral 
atom having only one optical electron with orbital angular mo¬ 
mentum 0 < / < 3, i.e. s-, p—, d- and /-states. We note 
that while the theories can be used for singly ionised atoms, 
this requires the calculation of the parameter E p for each state. 
For neutrals, this parameter can be approximated by a constant 
value in all cases, E p = -4/9 atomic units (see, e.g., Anstee 
1992, Barklem & O’Mara 1998). For this reason, the current 
work is restricted to neutral atoms. In the ABO broadening the¬ 
ory, spin is ignored. However, the calculation of depolarisation 
requires that spin be accounted for and thus we must consider 
specific groups of the periodic table. We treat the case of the 
alkaline earth group where the total angular momentum J-l. In 
addition we study the case of the alkali metals where the spin 
s = 1/2 and J = l + s. Therefore, in this paper we consider the 
cases of the p-states where J = 1=1, J= 1/2 or 3/2 and 7-states 
where J = 1=2, 7=3/2 or 5/2 and /-states where J=l=3, 7=5/2 
or 7/2. The case of complex atoms and atoms with hyperfine 
stmcture is considered in this work and its treatment is easily 
derived from the results concerned with simple atoms. This will 


be explained in sections 4.4 and 4.5 


Our goal is to retrieve a numerical relationship between the 
collisional depolarization rates and the collisional broadening 
line widths; that is D k (J)/w and C k {J —» J')/w. If one assumes 
that the velocity exponents are the same in the case of both 
broadening and depolarization, i.e. Aj_* jt = A, which is approx¬ 
imately the case, one finds that: 


D k (J) _ o- k (v 0 ) 
w o- w (v 0 ) ’ 
C k (J -> J') = ^>o) 
w cr w (v 0 ) 


( 11 ) 


3. Optimization approach 

To obtain these relationships, multipart data analysis and artifi¬ 
cial intelligence (Al) techniques are needed. Al techniques are 
becoming very useful as alternate approaches to conventional 
ones (Whiteson & Whiteson 2009). Within Al, genetic pro¬ 
gramming (GP) is a global optimization algorithm and an auto¬ 
matic programming technique that has been applied in physics 
and astrophysics (Cohen et al. 2003, El-Bakry & Radi 2007, 
Teodorescu & Sherwood 2008, El-dahshan 2009, Schmidt & 
Lipson 2009, Indranil et al. 2013). GP is a powerful tool that 
can be used to solve complex fitting problems. It is a rather 
recently developed evolutionary computation method for de¬ 
riving analytical functions and data analysis. GP is based on 
Darwin’s theory of evolution and uses population of individu¬ 
als, selects them according to an accuracy criterion, and pro¬ 
duces genetic variation using one or more genetic operators 
(Koza 1992). 

GP stores the individuals which are then presented as an 
expression tree for evaluation (Oakley 1994). GP provides the 
possible solutions to a given optimization problem, using the 
Darwinian principle of survival of the fittest. It uses biolog¬ 
ically inspired operations like replication, recombination and 
mutation. In GP, the chromosome or genome composed of a 
linear, symbolic string of fixed length, each consists of one or 
more genes. Each gene can be represented by an algebraic ex¬ 
pression and are composed of a head and a tail. The head con¬ 
tains symbols that represent functions F (with F=sqr,+,-,/,*, or 
other operators) and terminals T (with T=x,y,z, or other ran¬ 
dom constants), while the tail contains only terminals. 

As an example, the algebraic expression: 
((x*(x*y))+(x*y)) can be represented as a tree (Figure 
1). Note that {+, -, *) represents the head, and {x, y, x, z) 
represents the tail. This is the straightforward reading of 
the tree from left to right and from top to bottom. The set 
of functions and set of terminals must satisfy the closure 
and sufficiency properties. The sufficiency property requires 
that the set of functions and the set of terminals be able to 
express the solution of problem. The function set may contain 
standard arithmetic operators, mathematical functions, logical 
operators, and domain-specific functions. The terminal set 
usually consists of feature variables and constants. 

To measure the error between the actual data found by ABO 
(denoted by Z) and the data predicted by the numerical model 
(denoted by z mo dei ), we calculated the root of square mean error 
(RS ME): 


RSME = ^ 


n 

2 (^i — Z,model,i)~ 
1=1 


n 


( 12 ) 


where n is the number of the values of the data. The best solu¬ 
tion, which is adopted to obtain the results of this paper, corre¬ 
sponds to the minimum value of RS ME (Medhat et al. 2002). 


In the next sections we use the tabulated cross sections by ABO 
and DSB, for a wide range of cases and effective quantum num¬ 
bers n*, to establish relationships between cr w (v'o), cr^(vo) and 
0*j=>A v o). 


4. Results and discussion 

ABO presented the broadening cross sections in tabular form 
for general s - p, p - d, d - f and f - d transitions. Each 
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Fig. 4. Tree representation of the equation +(*(x,y),*(x,*(x,y))) 
i.e. ((x*(x*y))+(x*y)). 

cross section was obtained for given effective principal quan¬ 
tum numbers associated to the upper and lower levels of the 
transition. For instance, for p - d transitions, broadening cross 
sections cr w are tabulated, for a relative velocity of 10 4 m s 1 , 
with effective principal quantum numbers n p and n* r One can fit 
this table to obtain an analytic function of two variables cr w (n* p , 
n* d ). 

DSB tabulated the variation of the depolarization cross sec¬ 
tion cr k associated to the p, d. and /-states, for a relative veloc¬ 
ity of 10 4 m s 1 , with effective principal quantum numbers n*, 
n * d , and For instance, for the rZ-state, it is possible to obtain 
an analytic function of one variable cr k (n* d ). 

The ABO and DSB results are combined in order to ob¬ 
tain analytic functions of two variables: cr w (n *, cr k (n* p )), cr w (n p , 
cr k (n* d )) and cr w (n* d , cr^n*)), for the cases p, d and / states, rep- 
sectively. Note that for a transition between a lower level (with 
effective quantum number n*, ) and an upper level (with effective 
quantum number n*), our choice in this paper was to consider 
cr w as a function of // and n* and possible depolarization and 
polarization transfer cross sections cr k associated to n*. It is 
also possible to consider cr k associated to n*, instead of to n* u . 
However, typically in models it is the broadening associated to 
(nl, n* u ) and depolarization associated to the upper level n* r that 
is of greatest interest. 

In the following we use GP to treat each possible case and 
give the corresponding analytic functions. We note that the 
most usual situation in astrophysically modelling is that one 
is able to obtain and estimate for the collisional broadening 
rate, but requires an estimate for depolarization and polariza¬ 
tion transfer rates. 

4.1. s-p and p-s transitions 

For transtions between 5-states (Z = 0) and p-states (Z = 1), 
the cross sections cr w are given in Table 1 of Anstee & O’Mara 
(1991) as functions of the effective quantum number n* (corre¬ 
sponding to the s-state) and of the effective quantum number n* 
(corresponding to the p-state). For Z = 1, the non-zero depolar¬ 
ization and polarization transfer cross sections are cr 2 J= t , cr 2 =} p 


and cr°_ i/ 2 -, 3 / 2 - These cross sections are given as functions of 
n* in Table 6.1 of Derouich (2004). Note that we consider only 
the tensorial orders k — 0 and k — 2 which are relevant to the 
second solar spectrum studies (i.e. the scattering linear polar¬ 
ization spectrum). We wish to establish the relationships be¬ 
tween the function cr w (n*,n * p ) and the the functions cr 2 =l (n p ), 
o-} =3 / 2 («p and cr 2 J=l/2 ^ V2 (n* p ). Since o-} =| («p, <x® =1/2 ^ 3/2 (n*) 
and <T 2 J=1 p^ 3 p(n p ) depend on n p , the problem becomes to de¬ 
termine the general functions <x vv (n*,<xj =| ), or„(n*,cr 0 J=l/n ^ 3 p), 
and cr M ,(n*,cr5 =1/2 _ 3/2 ). 

4.1.1. Depolarization rate: J = J'= 1 and k = 2 (5=0) 

In the case of the depolarization rate where the total angular 
momentum J-l- 1, the relationship obtained from our GP mod¬ 
eling is: 

Z = -7. + 4. x X/5. + 35./Y 2 + (245./X) x (35. - Y) + 

(7 x 75. - X/7.) x (F + 5./X) x (Y - 2./Y), (13) 

where Z, 7, and X are introduced for the sake of clarity where, 
Z = cr w (n*,cr 2 =l ), Y = n*, andZ = cr 2 =] (n p ). 

To test the performance of our fitting method we employ 
“hare & hound” approaches consisting of comparison between 
the exact solution (called input) and the solution provided by 
the equation ( fTJfr (called output). Figure 2 represents, with a 
solid line, perfect correspondence (output=input) and the result 
of the fit are showed with squares. The averaged relative error 
is ~ 5%, which is sufficiently precise for applications. 



Fig. 5. To test the performance of our fitting method we employ 
“hare & hound” approaches consisting of comparison between 
the exact solution (called input) and the solution provided by 
equation ( fl3| (called output). The line shows the one-to-one 
relation. 

As an example, we apply Eq. ( fl3] > to determine the rela¬ 
tionship between the collisional depolarization and the colli¬ 
sional broadening rates associated to the photospheric Sr 1 4607 
A line. The 4607 A line is the resonance line of Sr i, namely 5 s 2 
'So —> 55 5 p l Pi, and for this line: n*(Srl)=2.13 and n*(Sri 
)=1.55. According to the DSB method (see Derouich 2004 and 
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pape rs therein), cr 2 =] =483 a.u. (atomic units). By applying Eq. 
(13 1 , we find that cr M ,(Sri )=430 implying that cr 2 =1 = 1.12 x cr w 
(r w = 1.12). We notice that this value is rather similar to the 
one adopted by Faurobert et al. (1995), where they assumed 
that cr 2 =1 = cr w (r w = 1) for the Sri 4607 A line. 


As mentioned, the more usual case is that one has access 
to w or (r w and wants to know D 2 or cr 2 =1 . An equation giving 
X ( X = cr 2 =1 ) as a function of Z ( Z-cr w ) and Y (Y — «*), 
was therefore determined using our GP model and the tables of 
collisional data given by ABO and DSB, namely: 


X = (2. x (Z/7. + 1 ,/Y) + Z-Y- 6.15) - 

(-2. xf + 2 x F 2 + 5.) x F 4 (14) 


In the case of the Sri line discussed above, one finds cr w (Sri 
)=405 a.u. by interpolation in the tables of Anstee & O’Mara 
(1995). Since Z=405 a.u. and T=1.55, using the general equa¬ 
tion (14 1 one can recover cr 2 =1 = 480.6 a.u., a value about 0.5% 
lower than the value calculated directly by DSB (see Derouich 
2004), 483 a.u. Thus the depolarization and polarization trans¬ 
fer rates can be obtained via the determination of the broad¬ 
ening rates. Note that by determining the depolarization and 
polarization transfer rates of simple atoms like Sri, the rates 
associated to complex atoms and atoms with hyperfine struc¬ 
ture can be easily obtained (see Sections 4.4 and 4.5 ). 

In the same homologous group, one can apply our results 
to the corresponding transitions in the Mgi and Cai atoms. 
For Mg i 3s 2 and 3 p 1 P\ levels one has n*(Mgi)= 1.33 and 
n * p (Mg i)=2.03. Then, using Derouich (2004), one has <x 2 =1 
=422 a.u. at n* p = 2.03. The application of Eq. (10) permits us 
to determine cr w (Mgi )=361 a.u. and thus <x 2 =1 (Mgi)=1.17 x 
cr w (Mgi) ( r w = 1.17). Similarly, for the case of Cai where 
n*(Cai )=2.07 and «*(Cai)=1.50, one has cr 2 =1 =444 a.u. and 
thus Eq (13 i gives cr w =394 a.u.. Therefore, cr 2 =] = 1.13 x cr w 
(r w = 1.13). I he results obtained for the resonance lines of Sri 
, Ca i and Mg i, are summarized in Table 1. 


4.1.2. Depolarization rate: J = J'=3I2 , k = 2 {s=M2) 

We adopt the same strategy as for J = J'-l in order to deter¬ 
mine the ratio r„ = D u ~ 3 t 2) - J =v 2< - l0> obtain the follow- 

W W OV(v o) 

ing relationship: 

Z = Y 4 x (2. + Y) + 49./Y 2 + 7. x X/(Y + 5.) + 
(-588. x Y + 98. x 7 3 + 294. x Y 2 ) x (5. + Y) x (7./X) (15) 

where Z = cr w (n*, cr 2 _ 3 ^ 2 ), Y = n*, and X = <r 2 J=3 / 2 (n* p ). For the 
Na I 3 s 2 S 1/2 and 3 p 2 Pi /2 levels one has «*(Na I)=1.63 
and n*(Na I)=2.12. According to Derouich (2004), cr 2 =3( , 0 =337 
a.u.. Thus, by applying Eq. ( fl5) , one finds that cr w = 433 a.u. 

Consequently, r w = Dk(J ^ ,2) = = 0.777 < 1. In 

addition, according to our GP method, the equation giving 
X = cr 2 J=i / 2 as a function of Z = cr w and Y = ri* is: 

X - Z - (Z/Y + 5.)/(5. x f) - 5. x f 5 (16) 


4.1.3. Transfer of population rate, 7=1/2, J'=3!2, k = 0 


The transfer rate of population is denoted by C k (J —> J') and 
the cross section of the transfer of population is . The 
relationship between cr w and CTj^j, is given by: 

Z = (10. - Y 2 )/(Y + 5.) x ((X/5. + 1. - Y) + (X - 56.)) 

+(1715. x Y 3 ) x l./X + (Y + 3.) x (7 3 ) x (Y + 9./4.) (17) 

where, Z = cr w , Y - «*, and X - (r k j _. t j,(n* p ). For example, the 
Na I 3 s 2 S 1/2 and 3 p 2 Pt ,/2 levels one has n*(Na I)=1.63 
and «*(Na I)=2.12 and using Derouich (2004) one finds that 

' _ 

thus r w - c U Y’ J 1 = ,Tj Y 1 ,' { '° ) = 0.6. Furthermore, by applying 


cr^ ,,=260 a.u.. The application of Eq. (1 17b gives cr w - 433, 


O\,.0o) 

qu; 

of Z = (T w and Y = n* is: 


our GP method, the equation giving X = cry =3 / 2 as a function 


X = Z/3. - Y 2 x Z/14. + 3. x Z/(3./T 2 + 5.) 


(18) 


4.2. p-d andd-p transitions 

4.2.1. Depolarization rate: J = J’= 2, k = 2 (s= 0) 

The depolarization rate for the cZ-states (1-2) have been calcu¬ 
lated by Derouich (2004) and references therein. In addition 
the collisional broadening by neutral hydrogen was obtained 
by Barklem & O’Mara (1997). After using our optimization 
method we obtain: 

Z = (7. x X)/(5. + Y) - (1./7.) x ( X/Y) - Y 3 x 3.5- 

490. x Y 2 /((X/5. - 49.) x (5. - Y) x (3,/5. - 2./Y)) (19) 

and 


X = (Z/7. + Z + 2.)/(30 ,/Y) + (Z + 4. x T) ■ 
x(Z + 5. - Y + 5./Z) 


■ (( 2 . 


Y)/ 12.) 
( 20 ) 


where, Z = cr w , Y = n*, and X = cr 2 =2 (n*). As an example, 
let us consider the case n* = 2 and n*.= 3. For n*.=3, Derouich 


(2004) found that o~ J=2 = 900 a.u. Equation (19 1 gives cr H ,=820 


a.u. The value given in Table 1 of Barklem & O’Mara (1997) 
is cr w =834 a.u., the relative error is 1.65 %. Figure 3 shows 
the comparison between the exact solution (the input from 
Barklem & O’Mara 1997) with solid line and the solution pro¬ 
vided by the equation ( p~9] > (called output) which is presented 
in squares. 


One should remark that the precision of equations (19 1 
and m are practically the same since the same optimization 
method and the same data are used to establish these two equa¬ 
tions. This remark is true for all equations of this paper giving 
Z and X. For this reason we usually check the accuracy of the 
equation giving Z. The reader could verify easily that the ac¬ 
curacy of the equation giving X is practically the same as the 
accuracy of the equation giving Z. 


4.2.2. Depolarization rate: J = J’=3I2, k = 2 
The relationship between cr w and cr 2 =3 ^(n*) is: 

Z = X x 2./(3. x Y) + 247. - 49. x Y + X+ 
(245. x Y x (7. - Y - T 2 ))/(-X/7. + Y x 6. + 16.) 


( 21 ) 
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r lv (Sr i) 

r w ( Cai) 

r w (Mgi) 

1.12 

1.13 

1.17 


Table 1 . r w = D { J 11 = {T i <ln> f or resonance lines of Sr i, Ca i, and Mg i. 

W W D-,„(v o) ° 



0 500 1000 1000 2000 2500 

Input: original data 


Fig. 6. To test the performance of our fitting method we employ 
“hare & hound” approaches consisting of comparison between 
the exact solution (called input) and the solution provided by 
the equation ( f]~9| ) (called output). The line shows the one-to- 
one relation. 


4.2.4. Transfer of population rate: 7=3/2, 7'=5/2, k = 0 
We find the function giving the relationship between cr M , and 
cr y=5/2("P : 

Z = A/(7. x F) + (1. - F) x 49. + A - 56. - (F - 7.) x (A/7.) 
+(5. x F 3 + F 5 x 5.)/[A/(10. x Y) - 8.)] (25) 

and 

X = 37. + Z - 6./F - 3. x (Z - 6 .)/(Y + 5.) - 2 x Y/3. (26) 

where, Z = cr w , Y = n p , and X = cr^ 2 _ )5 ^ 2 (np. At n* p = 2 and 
n* d = 3, we have from Derouich (2004) cr 2 =5 ^ =512 a.u. which 
gives cr w = 821 a.u.. The relative error in the calculation of the 
cr w from the relationship of Eq. < |25} is 1.6 %. 

4.2.5. Transfer of alignment rate: 7=3/2, 7'=5/2, k = 2 

Z = (X 2 /(5. + F))/((l. + Y) x 6.) + 49. 

+X/Y + (X + 3.) x 6. + (7./F 2 + 24.)/(—18. + A/10.) (27) 


and 

X = Z x Y/(Y + 1.) - (35. x Y) 

+(5. x(Z-l.)/7.)/(14. +2. xY) 

+(Z x Y/(Y + 1.) - 3. x F)/(8. x Y + 21.) (22) 

where, Z = cr„,, Y = n p , and X = cr” =3 ^ 2 (n*). Taking again the 
case of n* p = 2 and n* d =3, we have cr 2 =3 ^ = 496 a.u. which 
gives cr w -199 a.u.. The relative error in the calculation of the 
cr w from the relationship of Eq. pTj ) is 4.2 %. 

4.2.3. Depolarization rate: 7 = 7'=5/2, k = 2 

Using the optimization method, we find that the relationship 
between cr M , and <x 2 =5 ^(«*) is: 

Z = A x 2./(3. x Y) + 247. - 49. x Y + X+ 

(245. x Y x (7. - Y - F 2 ))/(-A/7. + Y x 6. + 16.) (23) 

and 


and 


A = (Z - 3.)/(68./9. + 3 ,/Y) - 2./Y+ 10./7. x Y + 10. (28) 


where, Z = cr„, Y - n* p , and A = cr^ 2 _^ 2 («p. At n* p = 2 


and n* d - 3, according to Derouich (2004), c r 3 / 2 ~> 5 /' , ( w P = 
a.u. which gives cr w - 826 a.u.. The relative error in the calcu¬ 
lation of the cr H , from the relationship of Eq. <[27} is 1.0%. The 

rat ior w =^ = ^ = 0 . 12 5. 

w w CT w (vo) 


4.3. d-f and f-d transitions 

4.3.1. Depolarization rate: 7 = 7'=3, k - 2 (s - 0) 

The depolarization rate for the /-states (1-3) have been cal¬ 
culated by Derouich (2004) and references therein. The col¬ 
lisional broadening by neutral hydrogen was obtained by 
Barklem & O’Mara (1998). The analytical relationship be¬ 
tween cr w and cr 2 =3 («*) is: 

Z = (15. x F 4 )/(A + 2.) x (5. x Y + 35.) x (F 2 - (7. + F)) (29) 
+(F - 2.)/6. x (1. + A) x (3. — F) + 

10./9. xA+11,/9.-F 2 


A = Z + 33. - Z/(5. + F) - Z/(F + 3.) + (F + 7.) x F (24) and 


where, Z = cr w , Y - n* p , and A = o-j =5 ^ 2 (n* d ). At n* = 2 and 
n* d - 3, we have from Derouich (2004) cr 2 =5 ^ 2 = 598 a.u. which 
gives cr w = 827 a.u.. The relative error in the calculation of the 
cr„, from the relationship of Eq. (23 i is 0.89 % and r w = = 

o*j(y o) 


O"».(vo) 


= 0.72. 


A = Z- l./F + 8 x F + 41. 

-(Z-51.-3. xF)/(Fx(F+1./3.)) (30) 

where, Z = cr w , F = /;*, and A = cr 2 (n*). This equation gives 
the relationship between cr w and cr 2 for any line, i.e. for all 
typical values of n* d and n*. Let us take values of the effective 
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Input: original data 


Fig. 7. To test the performance of our fitting method we employ 
“hare & hound” approaches consisting of comparison between 
the exact solution (called input) and the solution provided by 
the equation (291 (called output). The line shows the one-to- 
one relation. 


quantum number, n * d -3 and n* = 4. The cr w given by Barklem 
& O’Mara (1998) is 1419 a.u.. For the case n* = 4., Derouich 


(2004) found that cr 2 (> l 'f- 4) =1365 a.u.. According to the Eq. 


one shows easily that cr w =1464 a.u.. The relative error 
in the calculation of the cross section via Eq. ( [29] ) is -3.1%. 
Figure 4 shows the comparison between the exact solution (the 
input from Barklem & O’Mara 1998) with solid line and the 
solution provided by the equation (291 (called output) which is 
presented in squares. 


It is worth commenting on the noticeable difference in scat¬ 
ter in Figs. 5, 6 and 7, and the somewhat surprising result that 
the smallest scatter is seen in Fig. 6, for p-d and d-p transi¬ 
tions, rather than the scatter say correlating with orbital angular 
momentum quantum number. The scatter is largely traceable 
to non-smooth behaviour in the input data, especially the tabu¬ 
lated broadening cross sections. Figs. 5, 6 and 7, correspond to 
the cases of s-p (and p-s), s-p (and p-s), and d-f (and f-d) 
transitions, respectively, and the behaviour of the cross sections 
with n* is displayed graphically in the relevant ABO papers. 
There we see that non-smooth behaviour is especially seen at 
the extremities (i.e. large n*) of the tabulated data. This occurs 
due to the larger impact of oscillations in the electronic wave- 
functions on the calculations. That is, the wavefunctions have 
approximately n* —l — 1 nodes, and so at around n* > 1 + 2 such 
effects start to be seen in the broadening cross sections. We 
note that in the case of s-p and p-s transitions, the results are 
tabulated for .v-states up to n* = 3 = / + 3. However, in all other 
cases of the ABO data, the tabulations stop at n* = l + 2. Thus, 
the rather unexpected behaviour of the scatter is attributed to 
the larger range of n* calculated in the s-p and p-s case. If 
we restricted the s-p and p-s case to ,y-states with n* = 2, we 
could reasonably expect significantly reduced scatter in Fig. 5; 
see Fig. 1, upper panel, of Anstee & O’Mara (1995), showing 
the behaviour is very smooth in this regime. 


4.3.2. Depolarization rate: J = 7'=5/2, k = 2 (/ = 3) 

The relationship between cr w and <T? /2 0r‘.) is: 

Z = 2.xX-Yxl.-(5.-Y)xYx5. 

~(Y - 7./Y) x (Y - 2.) x (A/7.) (31) 

+(F - 2.) 2 x (35. IX) x Y 6 


and 


X = (27./10. + 3. x Y) x (3 x 7/2. - 2.) + 
3./7. x Z + (Z + 7) x 7/35. 


(32) 


where, Z = cr w , 7 = n* d , and X - cr^ 2 (n^)-l59 a.u.. For n * d =3 
and «*= 4, we find that cr w =1428 a.u. and that the error in the 


calculation of the cross section via Eq. (31 1 is -0.65 %. 


4.3.3. Depolarization rate: J = J'= 7/2, k = 2 (l = 3) 

The relationship between cr w and CTy 2 (n*f) is: 

Z - X - (10. x Y 2 ) + (7. x X)/(2. x 7 + 2.) - (7 + 1.) 

x(7. xX)x (3. - 7)/((5. - 7) x (A/5. - 86.)) + 2. (33) 


For n* d - 3 and /;'/= 4, the cr K inferred from Eq. (33 i is 1457 a.u. 
implying that the percentage of the relative error is -2.7 %. 
Since cr* 4)=824 a.u., the ratio r w = ^-21 = °~' (ln) _ 

7=7/2 v J w w cr w (vo) 

0.565. In addition, we find that: 


X = (Z/2. + 4. + 2 x 7) - (Z + 5.)/(4 x 7 2 ) - (-2./Z + 2./Z 2 ) 


x5. x (Z + 3.) x (5. + 7 + 2. x 7) 


(34) 


4.3.4. Transfer of population rate: 7=5/2, 7'=7/2, k = 0 
(1 = 3) 

We find that: 

Z = 2401. x (7 9 )/(A x (X - 7.)) + (2 x A + 12.)/(7 + 1.) 

+A-9.x 7-23. (35) 


and: 

A = (3. x Z/2. + Z x 7 + 8.)/(14 x 7 + 82.) + (10./7)/(Z + 1.) 
+5 x 7 x (5. + 7) + (3. + Z)/2. (36) 

where, Z = cr M ,, 7 = n* d , and A = c r ^ 2 _ > 7 / 2 ( w /)- For n j -3 and 
n*.= 4, the <x H . inferred from Eq. d35j) is 1447 a.u.. Therefore, by 
compring this value to the actual value calculated by Barklem 
& O’Mara (1998), one finds the percentage of the relative er¬ 
ror is -2 %. We notice that, according to DSB, = 

4.)=962 a.u. < cr w . 

4.3.5. Transfer of alignment rate: 7=5/2, 7'=7/2, k = 2 
(1 = 3) 

Our final relationship is: 

Z = A x (A - 5.)/7. + (7. + A)/(7 - 2.) 

+(35.-A/3.) x ((35.-A/2.) 

+(35. - A/7) x (7. - A) x (2. - 7)/(70. - 10. x 7)) 


(37) 
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where, Z = cr w , Y = n* d , and X - cr^ /? _ <7p (»* f ). For n* d - 3. and 
4., the cr w inferred from Eq. (|37ji is 1476.6 a.u. implying 
that the percentage of the relative error is -4 %. cr*(«* = 4) = 
101. < cr w . In addition, we obtain: 

X = 7. x (8. - Y/Z) + Z/35. - 245.245/(Z/21. - 30.) (38) 

4.4. Case of atoms with hyperfine structure 

Numerical results presented in the previous subsections are 
concerned with simple neutral atoms for which the electronic 
configurations have only one valence electron, and no hyper¬ 
fine structure. The problem of atoms with hyperfine structure 
can be easily dealt with by using the fact that, in the frame¬ 
work of the frozen nuclear spin approximation, the hyperfine 
depolarization and polarization transfer rates are given as a lin¬ 
ear combination of the fine rates D k (J ) and C k (J —> J') (e.g. 
Nienhuis 1976 and Omont 1977). Note that, in typical solar 
conditions, the frozen nuclear spin approximation is well sat¬ 
isfied. The hyperfine splitting is usually much lower than the 
inverse of the typical time duration of a collision and therefore 
one can assume that the nuclear spin is conserved (frozen) dur¬ 
ing the collision. 

4.5. Case of complex atoms 

The electronic configuration of a complex atom (Fe, Ti, etc.) 
often has one valence electron outside an open subshell with 
non-zero angular momentum. We assume that only the va¬ 
lence electron is affected by collisions with hydrogen atoms 
(the frozen core approximation). For more details about our 
approach concerning complex atoms we refer to Section 2. 
of Derouich et al. (2005a). In that paper and in papers that 
followed (Derouich et al. 2005b, Derouich & Barklem 2007, 
Sahal-Brechot et al. 2007), it has been shown that the expres¬ 
sion for the depolarization rate of a complex atom can be writ¬ 
ten as a linear combination of the depolarization rates of sim¬ 
ple atom. Thus, the relationships obtained and illustrated in this 
work are applicable to the case of complex atoms provided that 
one takes into account the formulae giving the depolarization 
rates of complex atoms as a function of depolarization rates of 
simple atoms. 

5. Conclusion 

A coherent scattering theory which accounts for partial or com¬ 
plete frequency redistribution in the presence of isotropic colli¬ 
sions is a challenging problem. Consequently, models of scat¬ 
tering are usually developed in a collisionless regime (e.g. 
Casini et al. 2014 and references therein). Stenflo (1994) pro¬ 
posed an approximate model including redistribution effects 
and accounts self-consistently for collisions. Chapter 10 of 
Stenflo (1994) confronted the problem of including properly 
the collisional depolarization rates and the collisional broaden¬ 
ing coefficients. 

In this paper, using GP techniques, we obtained accurate 
relationships for the collisional broadening as functions of the 


ddepolarization and polarization tranfer rates, accurate to typ¬ 
ically of order 5%. Since the the more usual case is that one 
has access to the collisional broadening and wants to know the 
depolarization and polarization tranfer rates, we determined ac¬ 
curate relationships giving the depolarization and polarization 
tranfer cross-sections as functions of the broadening line width. 
As such, these relationships should be useful to the solar polar¬ 
ization community. 

Interestingly, having the values of broadening coefficients 
and our relationships it is possible find the depolarization and 
polarization transfer rates for any hyperfine or fine-structure 
level of simple and complex atoms. 
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